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Genetic Analysis Systems and Methods 

BACKGROUND OF THE INVENTION 

The DNA that makes up human chromosomes provides the instructions that direct 
the production of all proteins in the body. These proteins carry out the vital functions of 
life. Variations in DNA often produce variations in the proteins, thus affecting the 
function of cells. Although environment often plays a significant role, variations or 
mutations in DNA are directly related to almost all human diseases, including infectious 
disease, cancer, inherited disorders, and autoimmune disorders. Moreover, knowledge of 
human genetics has led to the realization that many diseases result from either complex 
interactions of several genes or from any number of mutations within one gene. For 
example, Type I and II diabetes have been linked to multiple genes, each with its own 
pattern of mutations. In contrast, cystic fibrosis can be caused by any one of over 300 
different mutations in a single gene. 

Additionally, knowledge of human genetics has led to a limited understanding of 
differences between individuals when it comes to drug response — the field of 
pharmacogenetics. Since the first correlation over half a century ago of adverse drug 
responses with amino acid variations in two drug-metabolizing enzymes, plasma 
cholinesterase and glucose-6-phosphate dehydrogenase, careful genetic analyses have 
linked sequence polymorphisms in over 35 drug metabolism enzymes, 25 drug targets 
and 5 drug transporters with compromised levels of drug efficacy or safety. In the clinic, 
such information is being used to prevent drug toxicities; for example, patients are often 
screened for genetic differences in the thiopurine methyltransferase gene that cause 
decreased metabolism of 6-mercaptopurine or azathiopurine. Yet only a small percentage 
of observed drug toxicities have been explained adequately by the set of pharmacogenetic 
markers validated to date. In addition, insufficient therapeutic efficacy or unanticipated 
side effects in "outlier"" individuals when administered drugs previously demonstrated to 
be both safe and efficacious in clinical trials is a tremendous problem for health care 
practitioners and presents a significant dilemma to the pharmaceutical industry. 
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Disease-related and pharmacogenetic gene validation relies on elements of 
population and quantitative genetics and robust statistical metrics; however, the first step 
normally relies upon identification of a candidate target gene. To date, various 
biotechnological methods have been employed to identify candidate genes. For example, 
5 differential gene expression has been employed, essentially looking for differences in 
gene expression between affected and unaffected individuals or between treated and 
untreated individuals. In addition, protein-protein interaction maps to identify drug 
receptors and their immediate effectors have been used. Another approach involves 

" n mining human sequence databases for sequences similar to accepted disease-related or 

y 10 pharmacokinetic or pharmacodynamic regulators. 

4= 

fU Because any two humans are 99.9% similar in their genetic makeup, most of the 

%2 sequence of the DNA of their genomes is identical. However, there are variations in 

® DNA sequence between individuals. For example, there are deletions of many-base 

O stretches of DNA, insertion of stretches of DNA, variations in the number of repetitive 

H 15 DNA elements in noncoding regions, and changes in single nitrogenous base positions in 
the genome called "single nucleotide polymorphisms" or "SNPs". It is estimated that 
there are 3 to 4 million common SNPs that occur in at least 10 percent of people. These 
common SNPs do not occur independently but are passed from generation to generation 
in variable-length blocks of multiple SNPs, forming patterns across the genome. Such 
20 blocks of SNPs are called SNP haplotype blocks herein. 

The candidate gene identification strategy most relevant to SNPs is whole-genome 
association on various populations of individuals — that is, scanning the entire genomes of 
populations of individuals to correlate SNPs to disease or drug response. Such whole- 
genome analyses would provide a fine degree of genetic mapping and pinpoint specific 
25 regions of linkage. Methods have been proposed and are used in connection with whole 
genome analysis. For example, the methods described in U.S. Ser. No. 60/327,006, filed 
October 5, 2001, "Identifying Human SNP Haplotypes, Informative SNPs and Uses 
Thereof," assigned to the assignee of the present invention (Attorney docket no. 054801- 
5001-P2; 1005-P3, incorporated herein by reference for all purposes) have been proposed 



3 



Attorney Docket No. 1016 



for use in such applications. U.S. Ser. No. unassigned, filed September 18, 2001, 
"Human Genomic Polymorphisms", assigned to the assignee of the present invention 
(Attorney Docket no. 101083 19/1013P-1) incorporated herein by reference for all 
purposes, provides the identity of SNPs and SNP haplotype blocks across one 
representative chromosome, ie. Chromosome 21. 

While meeting with success, it is desirable to increase the speed and efficiency at 
which such analyses can be performed, as well as to decrease the cost of performing such 
analyses. 

SUMMARY OF THE INVENTION 

Improved methods, devices, and systems for performing genetic analysis are 
provided. According to one embodiment of the invention a method is provided in which 
genetic DNA or derivatives therefrom is scanned from a plurality of individuals. The 
scanning uses, for example, a gel based sequencer or high density array to identify genetic 
polymorphisms in a plurality of individuals. For the polymorphisms that occur in a 
plurality of individuals, but not all of the scanned bases, genetic DNA or derivatives 
therefrom from additional individuals is evaluated to identify which of the variants occur 
in the additional individuals. Based the results of scanning the additional individuals, the 
method provides for the identification of blocks of the polymorphisms that tend to be 
inherited together. The blocks of polymorphisms can then be used in a wide variety of 
applications such as drug target identification, identification of diagnostics, clinical trial 
evaluation, pharmacogenomics, and other applications. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

Throughout the disclosure various patents, patent applications and publications 
are referenced. Unless otherwise indicated, each is incorporated by reference in its 
entirety for all purposes. 
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U.S. Ser. No. 60/327,006, filed October 5, 2001, "Identifying Human SNP 
Haplotypes, Informative SNPs and Use Thereof," assigned to the assignee of the present 
invention (Attorney docket no. 054801 -500 1-P2) incorporated herein by reference for all 
purposes, describes methods for identifying single nucleotide polymorphisms, assigning 
such polymorphisms to groups, and using such groups in the drug discovery process. Fig. 
1 provides an overall diagram illustrating such processes. In Fig. 1, variant locations in a 
genome, such as the human genome, are identified in step 110. The processes used in 
step 110 can include, for example, scanning genomic regions for variants using 
DNA/RNA arrays such as those available from Affymetrix, Inc. (such as discussed in 
general in U.S. Pat. No. 5,861,242, incorporated herein by reference for all purposes), the 
use of conventional gel sequencing techniques such as those based on gel sequencers 
available from Applied Biosystems, Inc., capillary sequencing techniques, or other 
techniques. Preferably, only common SNPs, such as those that occur in more than 10% 
of the population in question, are used in the analysis since, for example, these SNPs will 
represent evolutionary older variants and, therefore, tend to correlate with more 
common diseases/states in a larger fraction of the population. 

In preferred embodiments, the genomic DNA (or a derivative therefrom, such as 
PCR DNAs or RNA amplicons or the like) is scanned for variants using a high density 
DNA array. Using one technology, an array of probes is formed on a high density array. 
The array of probes includes a) probes that are complementary to portions of a baseline 
sequence of the genetic DNA (such as that obtained from a first pass sequencing effort) 
and b) one or more probes that are complementary to variants of the baseline sequence of 
the genetic DNA. The genetic DNA, or a derivative thereof, is then labeled and 
hybridized to the array and scanned. Based on the locations where hybridization is 
detected one can detect whether the bases in the genetic DNA are "wild type" or variant. 
In the cases where the particular base is the same as the reference sequence, the labeled 
material will hybridize more strongly to the wild type probe. In those cases where the 
DNA sequence varies from the reference sequence, the labeled material will hybridize 
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more strongly to one of the variant probes. Based on this information, one can call a 
particular base as being the same or variant from the reference sequence. 

In preferred embodiments, more than 10,000 bases are scanned for variants. In 
more preferred embodiments, more than lxlO 6 bases are scanned for variants, more 
5 preferably more than 2xl0 6 bases, more preferably lxlO 7 , and more preferably more than 
lxlO 8 , and more preferably more than lxlO 9 bases are scanned for variants. In preferred 
embodiments at least introns are scanned for variants, and in more preferred 
embodiments both introns and exons are scanned for variants. In preferred embodiments 
the scanned nucleic acids are genomic DNA. In most preferred embodiments, such DNA 
10 is from a mammalian organism such as a human. In preferred embodiments, more than 
10% of the genomic DNA from the organism is scanned, in more preferred embodiments 
more than 25% of the genomic DNA from the organism is scanned, in more preferred 
embodiments, more than 50% of the genomic DNA from the organism is scanned, and in 
O most preferred embodiments, more than 75% of the genomic DNA is scanned. 
J«j 15 In step 120 the variants are assigned to groups. Such groups of variants are blocks 

jf of variants (in most cases, single nucleotide polymorphisms or "SNPs") that tend to be 
inherited together. Such "blocks" or "SNP variant haplotype blocks" have likely arisen 
from an evolutionary point of view because the edges of blocks will tend to be the 
boundaries of ancient recombination events in the species. In step 130, representative 
20 variants from each block are used in an association study to identify associations of such 
blocks to a phenotype. Such association studies may be run using a variety of 
technologies such as a) Affymetrix genotyping arrays, b) "Invader" assays available from 
Third Wave, such as described in Agarwal, P. et al. "Comparison study for identifying 
promoter allelic polymorphism in interleukin 10 and tumor necrosis factor alpha genes", 
25 Diagn Mol Pathol 9, 158-64(2000); Cooksey, R.C., Holloway, B.P., Oldenburg, M.C., 
Listenbee, S. & Miller, C.W. "Evaluation of the Invader assay, a linear signal 
amplification method, for identification of mutations associated with resistance to 
rifampin and isoniazid in Mycobacterium tuberculosis" Antimicrob Agents Chemother 44, 
1296-301 (2000); and Mein, C.A. et al. "Evaluation of single nucleotide polymorphism 
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typing with Invader on PCR amplicons and its automation" Genome Res 10, 330-43 
(2000), all incorporated herein by reference, or c) 'Taqman" assays, such as described in 
Arnold, B.A., Hepler, R.W., and Keller, P.M. "One-Step Fluorescent Probe Product- 
Enhanced Reverse Transcriptase Assay" BioTechniques 25(1):98-106, 1998; Becker K., 
5 D. Pan and C.B. Whitley. 1999. Real-time quantitative polymerase chain reaction to 
assess gene transfer. Hum. Gene Ther. 70:2559-2566, 1999; Berg, T., Miller, A.R., Platz, 
K.P., Hohne, M., Bechstein, W.O., Hopf, IL, Wiedenmann, B., Neuhaus, P., and Schreier, 
E. "Dynamics of GB virus C viremia early after orthotopic liver transplantation indicates 
H extrahepatic tissues as the predominant site of GB virus C replication" Hepatology 
Q 10 29(l):245-249, 1999; or a variety of other technologies. When array based technologies 
nj are used, as an example, one makes arrays similar to those described above, except that 

only the variant positions are "tiled" onto the array. 
*0 Once it has been determined which regions of the genome are associated with, for 

O example, a disease state, drug screening assays are run to identify one or more molecules 
:S 15 that mediate the relevant portions of the genome in step 140 to identify therapeutics 150. 
5 In addition, the associations themselves may be used directly in genetic diagnostic assays, 
m immunoassays or other diagnostic tests in step 160. In alternative embodiments the 
variant regions are used in studies to evaluate drug response such as in pharmacogenomic 
studies. The results are then used to stratify patents according to drug response. The 
20 stratification may, in some embodiments, be used to stratify individuals that may exhibit a 
toxic response from those that will not. In other cases, the stratification may separate 
those that will be non-responders from those who will be responders. 

Fig. 2 provides one illustration of showing how variants, usually SNPs, occur in 
groups or blocks in a genome. Such SNP haplotype blocks are chromosomal regions that 
25 tend to be inherited as a unit, with a relatively small number of common forms. Each line 
in Fig. 2 represents portions of the haploid genome sequence of different individuals. As 
shown therein, individual W has an "A" at position 241, a "G" at position 242, and an 
"A" at position 243. Individual X has the same bases at positions 241, 242, and 243. 
Conversely, individual Y has a T at positions 241 and 243, but an A at position 242. 
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Individual Z has the same bases as individual Y at positions 241, 242, and 243. The 
variants will most commonly be biallelic, ie. occur in only one of 2 forms, not all four 
possibilities. Variants in block 261 will tend to occur together. Similarly, the variants in 
blocks 262 and 263 will tend to occur together. Of course, only a few bases in a genome 
5 are shown in Fig. 2. In fact, most bases will be like those at position 245 and 248, and 
will not vary from individual to individual. 

This tendency of SNPs to occur together in SNP haplotype blocks allows 
for a single or a few SNPs in a block to be used in an association study, rather than all of 
the SNPs in that block. For example, by evaluating only the SNP at position 241, the 
10 SNPs at positions 242 and 243 can be predicted without performing an assay to identify 
SNPs 242 and 243. If position 241 contains an A, position 242 will contain a G and 
position 243 will contain an A. Conversely, if position 241 contains a T, positions 242 
and 243 respectively will contain an A and a T respectively. Similarly, a SNP in block 
262 can be used to predict the other SNPs in block 262. For example, if SNP 244 is a T, 
15 it will be known that position 246 is a C, and position 247 is a G. Conversely, if position 
244 is a C, position 246 is an A, and position 247 is a C. 

According to the invention herein, all or substantially all of the bases in a 
genomic region are scanned for common SNPs in a first, limited, sample set. Then, using 
substantially only those SNPs identified in the first scanning step, additional individuals 
20 are genotyped at those SNP locations that are identified from a full scan of the genomic 
region of interest. The first sample set is sufficiently large that the full scan of genomic 
DNA identifies a sufficient number of SNPs with a sufficiently high level of confidence, 
but is not large enough to necessarily identify all of the block boundaries to a necessary 
level of confidence for the particular experiment. The second sample set is of a sufficient 
25 size that it will allow definition of block boundaries with a sufficient level of confidence 
for the experiment. 

Fig. 3 illustrates the process in greater detail. As shown, at step 310 a subset of 
"n" genomic DNAs are scanned across all or a portion such genomes at every base 
location for common SNPs or variants. Thereafter, using genotyping technologies such 
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as those described above, only those variant locations identified in step 310 are used to 
genotype an additional "m" individuals. Using the genotypes accumulated from steps 
310 and 320, variant blocks are assigned. Using the invention herein, the cost and burden 
of performing genetic analysis, especially whole genome analysis, is significantly reduced 
since, for a fraction of the population analyzed, only a fraction of the bases of the 
genomic DNA of interest needs to be evaluated in assigning variants to respective blocks. 

In preferred embodiments, copies of more than 10 different monoploid genomic 
DNAs are scanned for variants at step 310 to be used in the further genotyping analysis 
(i.e. "n" is 10 or more). In more preferred embodiments more than 15 different copies of 
the genomic DNA are scanned for variants to be used in the further genotyping analysis 
(i.e. "n" is 15 or more). In more preferred embodiments more than 18 copies of the 
genomic DNA are scanned for variants to be used in the further genotyping analysis (i.e. 
"n" is 18 or more). In more preferred embodiments more than 20 genomic DNAs are 
scanned for variants to be used in the further genotyping analysis (i.e. "n" is 20 or more). 
In more preferred embodiments more than 25 copies of the genomic DNA are scanned for 
variants to be used in the further genotyping analysis (i.e. "n" is 25 or more). 

In preferred embodiments, the full scanning of genomic DNA is followed up by 
genotyping at step 320 of additional samples to evaluate representative SNPs discovered 
in the full scans. In preferred embodiments, the genotyping step 320 analyzes more than 
10 different monoploid genomic DNAs (ie. "m" is 10 or more). In more preferred 
embodiments, the genotyping step 320 analyzes more than 15 different genomic DNAs. 
In more preferred embodiments, the genotyping step 320 analyzes more than 20 different 
genomic DNAs. In more preferred embodiments, the genotyping step 320 analyzes more 
than 25 different genomic DNAs. In more preferred embodiments, the genotyping step 
320 analyzes more than 30 different genomic DNAs. 

Fig. 4 illustrates the outcome of the process of Fig. 3 in a hypothetical situation. 
The top portion of Fig. 4 illustrates the sequence of a hypothetical stretch of DNA, with 
the variant positions indicated and variant block boundaries drawn. Of course, the top 
portion of Fig. 4 illustrates the stretch of genomic DNA with mutations marked and the 



Attorney Docket No. 1016 



block boundaries drawn, although this information would not be known ab initio. 
Sequencing results 410 show the results of sequencing haploid DNA of 3 individuals. As 
shown, in general it is possible to have identified a large fraction of the common SNPs 
after a relatively small number of individuals have been sequenced. In the case in Fig. 4, 
the SNPs at each location shown in the top portion of Fig. 4 have been identified, as 
indicated by check marks. However, if further individuals are not evaluated, the block 
boundaries would not be correctly identified at this stage. For example, while one could 
at this stage draw block boundaries between blocks 420 and 430 (note that the first C ->G 
variant predicts the first G->A variant, and the first C->T variant predicts the second 
C->T variant), it is not possible to distinguish between the blocks 430 and 440 at this 
stage. At this stage it appears that the first C->T variant would predict the first and 
second T-> A variants. Accordingly, a more statistically significant sample set is required 
to draw the block boundaries. However, since relevant SNPs have been identified, it is 
now possible to genotype only the variant positions to complete the process of identifying 
block boundaries without sequencing the entire stretch of genomic DNA. 

The results of performing a genotyping process on only the SNPs in another 
hypothetical genomic sample are shown in Fig. 4 at 420. As shown, by performing this 
additional genotyping step, it is now possible to see that blocks 430 and 440 are 
distinguishable. Specifically, it is now possible to see that the first C->T variant does not 
track with the first and second T->A variants, but instead, the first C->T variant can be 
used to predict only the second C->T variant and the first T->A variant can be used only 
to predict the second T-> A variant. 

Without wishing to be bound by any particular theory of operation, it is believed 
that the basis on which the assays may be successfully utilized in this manner is found in 
the fact that a few individuals may be used to find the most common SNPs. The most 
common haplotypes account for most individuals, but less common haplotypes could be 
mis-ordered in terms of frequency of occurrence if a larger sample set is not used. 
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Example(s) 

Using Affymetrix whole wafers, a region of genomic DNA from chromosome 21 
was scanned for polymorphic regions, generally using the methods of U.S. Ser. No. 
60/327,006, filed October 5, 2001, "Identifying Human SNP Haplotypes, Informative 
SNPs and Use Thereof," assigned to the assignee of the present invention (Attorney 
docket no. 054801-5001-P2; 10O5-P3, incorporated herein by reference for all purposes). 
19 additional copies of chromosome 21 were then evaluated using the same techniques. 
"Common" SNPs were evaluated, meaning those that occur with an allele frequency of 
more than 10% in this particular example. To evaluate whether one may use 19 scans for 
SNP detection, rather than 38, it was determined how many SNPs would produce an 
allele frequency of >10% across all 38 scans, but are not considered common 
polymorphism in a set of just 19 scans. To have an allele frequency of >10% out of 38, a 
SNP would need to have a frequency of >20% in 19 scans, if it is not polymorphic at all 
in the other 19 scans. In the analysis performed, 89% of common SNPs were identified 
by evaluating just 19 scans. SNPs identified in the above analysis are shown in Tables 1, 
2, and 3. Table 1 shows the SNPs that would have been found in both of the first and last 
19 scans of a 38 scan set. Table 2 shows the SNPs that would have been identified by 
evaluating the first 19 of the 38 scans of chromosome 21, but would not have been 
identified by evaluating only the second 19 scans. Table 3 shows the SNPs that would 
have been identified by evaluating the last 19 of 38 scans of chromosome 21, but would 
not have been identified by evaluating only the first 19 scans. There are 3170 SNPs 
identified in both the first and last 19 scans. There are 361 SNPs found in only the first 
19 scans, and 516 in only the second 19 scans. Therefore, there are ((0.5) * 
(361+516))/(3170+361+516) or 11% of the SNPs, on average, that would not be found in 
only half of the number of scans in chromosome 21. 
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Table 1 
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It is to be understood that the above description is illustrative and not restrictive. 
The scope of the invention should, therefore, be determined not with respect to the above 
description, but instead with respect to the appended claims, along with the full scope of 
equivalents to which the claims are entitled. 
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